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Abstract. The LIGO Scientific Collaboration recently reported a new upper limit on an 
isotropic stochastic background of gravitational waves obtained based on the data from the 3rd 
LIGO science Run (S3). Now I present a new method for obtaining directional upper limits 
that the LIGO Scientific Collaboration intends to use for future LIGO science runs and that 
essentially implements a gravitational wave radiometer. 



1. Introduction 

The LIGO Scientific Collaboration analyzed the data from the 3rd LIGO science run (S3) looking 
for an isotropic background of gravitational waves. In Ref. J s it was reported that a 90%- 
confidence Bayesian upper limit on f2 gw (/) of 8.4 x 10 -4 was achieved. Q gw (f) is the energy 
density per logarithmic frequency interval associated with gravitational waves normalized by 
the critical energy density p c to close the universe. However it is possible that the dominant 
source of stochastic gravitational waves in the LIGO frequency band comes from an ensemble 
of astrophysical sources (e.g. [HIE])- If such an ensemble also turns out to be dominated by its 
strongest members the assumption of isotropy is no longer very good. Then one should look for 
anisotropics in the stochastic gravitational wave background. This was addressed in Ref. 
but they focused on getting optimal estimates for each spherical harmonic. If the stochastic 
gravitational wave background is indeed dominated by individual sources one can get a better 
signal-to-noise ratio by obtaining optimal filters for small patches of the sky. 

I present such a directional method that essentially implements a gravitational wave 
radiometer. The algorithm has been implemented and will be used to analyze future LIGO 
science runs starting with S4. 



2. Search for an isotropic background 

The data from the first 3 LIGO science runs was analyzed with a method described in detail 
in Ref. ^ |2] . The data streams from a pair of detectors were cross-correlated with a cross- 
correlation kernel Q chosen to be optimal for an assumed strain power spectrum H(f) = H(\f\) 
and angular distribution P(&) = 1 (isotropic distribution). Specifically, with Si(f) and 82(f) 
representing the Fourier transforms of the strain outputs of two detectors, this cross-correlation 
is computed in frequency domain segment by segment as: 



/oo roc 
df / df'8 T {f-f)Sl(f)Q(f)S 2 {f) 
-00 J — 00 



(1) 



where 5t is a finite-time approximation to the Dirac delta function. The optimal filter Q has 
the form: 

Q(/) - A A(/)P 2 (/) (2) 

where A is a normalization factor, Pi and P2 are the strain noise power spectra of the two 
detectors, H(f) is the strain power spectrum of the stochastic background being searched for 
(= S gw (f) in Ref. 0121) and the factor 7i so takes into account the cancellation of an isotropic 
omni-directional signal (-P(f2) = 1) at higher frequencies due to the detector separation. 7i so 
is called the overlap reduction function jS] and is given by (the normalization is such that 
7iso(/ = 0) = l for aligned and co-located detectors): 

7iso(/) = ^E I d& Ff(Q)F 2 A (n) (3) 
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where Ax = x 2 — %i is the detector separation vector and 

Ff(Q) = e^-iXfXt - Y?yf) (4) 

is the response of detector i to a zero frequency, unit amplitude, A = +, x polarized gravitational 
wave. 

The optimal filter Q is derived assuming that the intrinsic detector noise is Gaussian 
and stationary over the measurement time, uncorrelated between detectors, and uncorrelated 
with and much greater in power than the stochastic gravitational wave signal. Under these 
assumptions the expected variance, Oy, of the cross-correlation is dominated by the noise in 
the individual detectors, whereas the expected value of the cross-correlation Y depends on the 
stochastic background power spectrum: 
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Here the scalar product (•,•) is defined as (A, B) = p 00 A*(f)B{f)P l {f)P 2 {f)df and T is the 
duration of the measurement. 

To deal with the long-term non-stationarity of the detector noise the data set from a 
given interferometer pair is divided into equal-length intervals, and the cross-correlation Y and 
theoretical cry are calculated for each interval, yielding a set {l/jOy^} of such values, with / 
labeling the intervals. The interval length can be chosen such that the detector noise is relatively 
stationary over one interval. In Ref. [2j the interval length was chosen to be 60 sec. The 
cross-correlation values are combined to produce a final cross-correlation estimator, Y op t, that 
maximizes the signal-to- noise ratio, and has variance cr^ pt : 

Y op t = T,i<Ty?Yi/(T^ t , 0-~pt = E/Oyf • (6) 

Since the LIGO Hanford and Livingston sites are separated by 3000km the overlap reduction 
function for this pair has already dropped below 5% around each interferometer's sweet spot 
of 150 Hz. One unfortunate drawback of this analysis thus is the limited use it makes of 
the individual interferometer's most sensitive frequency region. Moreover, if the dominant 
gravitational wave background would be of astrophysical origin the assumption of an isotropic 
background is not well justified. If for example the signal is dominated by a few strong sources 
a directed search can achieve a better signal-to-noise ratio. 



3. Directional search: a gravitational wave radiometer 

A natural generalization of the method described above can be achieved by finding the optimal 
filter for an angular power distribution P(Cl). In this case Eq. \5)p generalizes to 



where 7^ is now just the integrand of 7i so , i.e. 

7« = ^E^'*^(^(fi) (8) 
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and H(f) is the strain power spectrum of an unpolarized point source, summed over both 
polarizations. Note that 7^ also becomes sidereal time dependent both through Ax and F^(U). 

Eq. [7| was used in Ref. [5,, as a starting point to derive optimal filters for each spherical 
harmonic. However if one wants to optimize the method for well located astrophysical sources 
it seems more natural to use a P(£l) that only covers a localized patch in the sky. Indeed, for 
most reasonable choices of H(f), the resulting maximal resolution of this method will be bigger 
than a several tens of square degrees such that most astrophysical sources won't be resolved. 
Therefore it makes sense to optimize the method for true point sources, i.e. P(&) = 5 2 (ti,£l'). 

With this choice of P(&) the optimal filter Q^, for the sky direction 0' becomes 



Qfr(/)=* y," p,,' ( 9 ) 



and the expected cross-correlation Ya, and its expected variance oh„ are 

4 ft , = (^)-(^) 2 -^fe,Qn') > ^> =r (°ft"ftf) (10) 
3.1. Integration over sidereal time 

Since the non-stationarity of the detector noise also needs to be dealt with, processing the data 
on a segment by segment basis still makes sense. However 7^ changes with sidereal time. By 
setting it to its mid-segment value one can get rid of the 1st order error, but a 2nd order error 
remains and is of the order 
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with / the typical frequency and d the detector separation. For T seg = 60 sec, / = 2 kHz and 
d = 3000 km this error is less than 1%. 

Thus the integration over sidereal time for each 0,' again reduces to the optimal combination 
of the set {17,057}^, given by Eq. |SJ The only difference to the isotropic P(Cl) = l case is that 
the optimal filter Q^, is different for each interval I and each sky direction Q,' . 

3.2. Numerical aspects 

To implement this method one thus has to calculate 
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for each sky direction Q' and each segment /. This can be calculated very efficiently by realizing 
that 7^ splits into a DC part 1/2Ea^(^)^2 1 (^) an d a phasor exp(i27r/fi • Ax/c). For both 
integrals the DC part can be taken out of the frequency integration, leaving all the directional 
information of the integrands in the phasor. Thus, with N the number of sky directions f2', 
instead having to calculate 2N integrations, it is sufficient to calculate one fast Fourier transform 

and one integral per segment, and read out the cross-correlation Y^, at the time shift r = Cl' •—. 

Since the fast Fourier transform of SI* S2 H / (P1P2) is sampled at /sample = 2/N yqu ist it is 
necessary to interpolate to get the cross-correlation Y^, at the time shift r. However, by choosing 
a high enough Nyquist frequency and zero-padding the unused bandwidth this interpolation error 
can be kept small while the overall efficiency is still improved. 



3.3. Comparison to the isotropic case 

It is interesting to look at the potential signal-to-noise ratio improvement of this directional 
method compared to the isotropic method if indeed all correlated signal would come from one 
point Cl', i.e. {SIS2) = l&H. The ratio between the two signal-to-noise ratios works out to 

SNRi„ _ ggVg _ [7. S o,7n,] 

with [A, B] = (AiH/(Pi t iP2 t i),BiH/(Pi t iP2 i i)) and i the index summing over sidereal time. 
This ratio is bounded between —1 and 1, i.e. the directional search not only performs better in 
this case but, for a point source at an unfortunate position, the isotropic search can even yield 
negative or zero correlation. 

It is also possible to recover the isotropic result as an integral over the sky. The definitions 
of 7i so and 7^ (Eq. |3]and|HJ) imply 

*EV2r a =ljdn Y**?- 2 , a°r 2 = {i) 2 j^j ^y 2 ( 14 ) 

Here cr? ^, is the covariance of the 2 sky directions Q, and Cl 1 . It is given by 

^"jtfitoQ*) (15) 
and describes the antenna lobe of the gravitational wave radiometer. 



3.4- Achievable sensitivity 

The l-cr sensitivity of this method is given by 

H a(/) = ^H(f) = (16) 

^V^y (f-oo 7 p 1 p 2 rf/)sidereal day 

H sens is somewhat dependent on the declination and in theory independent of right ascension. 
In practice though an uneven coverage of the sidereal day due to downtime and time-of-day 
dependent sensitivity will break this symmetry, leaving only an antipodal symmetry. 

For the initial LIGO Hanford 4km - Livingston 4km pair (Hl-Ll), both at design sensitivity, 
and a flat source power spectrum H{f) = const this works out to 

HST^^LSxIO-^Hz- 1 ^)* (17) 



with a 35% variation depending on the declination. This corresponds to a gravitational wave 
energy flux density of 



m 2 Hz V 100 Hz 

3. 5. Results from simulated data 

The described algorithm was implemented such that it is ready to run on real LIGO data. 
However, in order to test the code, the real data was blanked out and simulated Gaussian noise 
uncorrelated between the 2 detectors and with a power spectrum shape equal to the LIGO 
design sensitivity was added. Real lock segment start and stop times were used. This takes into 
account the non-uniform day coverage. To get a quicker turn-around during testing the code 
was only run on 1.7 days of integrated simulated data. The signal power spectrum was assumed 
to be flat, H(f) = const. 

The algorithm was run on a 360 x 181 point grid covering the whole sky. While this clearly 
over-samples the intrinsic resolution - for the H (/) = const case the antenna lobe has a FWHM 
area of 50 — 100 deg 2 , depending on the declination - it produces nicer pictures as final product 
as shown in Figure ^ 
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Figure 1. An example map of the signal-to- noise ratio SNR = l^ p t/<7 pt for simulated Gaussian 
noise (see text). The visible structure - fringes with opposite tilt on the northern and southern 
hemisphere as required by the antipodal symmetry of the antenna lobe - is due to the antenna 
lobe with which the random background is convolved. 

For Figure [21 additional coherent noise with a flat power spectrum H{f) and a sidereal time 
dependent time shift and amplitude modulation appropriate for a true point source was added. 
The data was produced by splicing together short pieces of data with fixed time shift. The 
injected source has a signal-to-noise ratio of 14 and is clearly recovered. This figure also shows 
the typical shape of the radiometer antenna lobe which is given by Eq. 1151 In particular it shows 
that areas adjoint to the main lobe get a negative correlation. This means that a particularly 
unfortunate set of sources could in principle cancel a lot of the signal. 
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Figure 2. A point source with signal-to-noise ratio (SNR) of 14 was injected at the position of 
the Virgo galaxy cluster (12. 5h, 12.7°). It is nicely recovered. This map also shows the typical 
structure of the antenna lobe including negative correlation in regions adjoint to the main peak. 



4. Conclusion 

The presented gravitational wave radiometer method aims for optimal detection of localized 
stochastic gravitational wave sources with a given power spectrum H (/) and significant uptime. 
It produces a map of point estimates and corresponding variances. The point estimate map 
corresponds to the true power distribution of gravitational waves P(£l) convolved with the 
antenna lobe of the radiometer and an uncertainty determined by the detector noise. This 
antenna lobe in turn depends on the assumed source power spectrum H(f). 

The method is well suited for searching for a stochastic gravitational wave background of 
astrophysical origin and the LIGO Scientific Collaboration intends to use it on future LIGO 
science runs starting with S4. 

The author gratefully acknowledges his colleagues in the LIGO Scientific Collaboration whose 
work on building and commissioning ground based gravitational wave detectors as well as 
computing infrastructure made this research possible. In addition he acknowledges the support 
of the United States National Science Foundation for the construction and operation of LIGO 
under Cooperative Agreement PHY-0107417. 

References 

[1] B. Abbott et al, "Upper Limits on a Stochastic Background of Gravitational Waves," Phys. Rev. Lett. 

accepted for publication 
[2] B. Abbott et al, Phys. Rev. D 69, 122004 (2004). 

[3] T. Regimbau, J. A. de Freitas Pacheco Astron. Astrophys. 376, 381 (2001) 

[4] L. Bildsten, Astrophys. J. L89-L93 (1998) 

[5] B. Allen and A.C. Ottewill, Phys. Rev. D 56, 545 (1997) 

[6] N. J. Cornish, Class. Quantum Grav. 18, 4277 (2001) 

[7] B. Allen and J.D. Romano, Phys. Rev. D 59, 102001 (1999). 

[8] N. Christensen, Phys. Rev. D 46, 5250 (1992); E.E. Flanagan, Phys. Rev. D 48, 389 (1993). 
[9] B. Barish and R. Weiss, Phys. Today 52, 44 (1999). 
[10] B. Abbott et al., Nucl. Instrum. Methods A 517/1-3, 154 (2004). 



